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Light hadron spectrum with two flavors of O(a) improved dynamical 
quarks : final results from JLQCD* 
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We present the final results of the JLQCD calculation of the light hadron spectrum and quark masses with 
two flavors of dynamical quarks using the plaquette gauge action and fully 0(a)-improved Wilson quark action 
at P = 5.2. We observe that sea quark effects lead to a closer agreement of the strange meson and baryon masses 
with experiment and a reduction of quark masses by 20-30 %. 



1. Introduction 

Precise determination of the light hadron spec- 
trum with dynamical up, down and strange 
quarks is a key step for validating QCD at low 
energies. As a step toward this goal, we have pur- 
sued a systematic calculation in two-flavor QCD 
using the non-perturbatively 0(a)-improved Wil- 
son quark action [jjj. In this article, we present 
the final results of the hadron spectrum and quark 
masses from this calculation. A test of the one- 
loop chiral perturbation theory formula in our lat- 
tice data is discussed in a separate talk || . 

Our simulations were performed at j3 = 5.2, 
where the scale is a = 0.0887(11) fin, on 12 3 x48, 
16 3 x 48 and 20 3 x 48 lattices. We simulated 
five sea quark masses i^sca = 0.1340, 0.1343, 
0.1346, 0.1350 and 0.1355, which cover the range 
of mps jS ea/ m v.sca = 0.80-0.60. For further de- 
tails of the simulations, we refer the reader to 
the previous reports Q and a forthcoming paper 
(3). The main progress after the last conference 
is that we accumulated doubled statistics (12000 
HMC trajectories) on the largest lattice. 
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2. Finite size effects 

Finite size effects (FSE) are more pronounced 
in full QCD than in quenched QCD. Since our 
largest volume size ~ (1.8 fm) 3 is not so large, it 
is important to check FSE in our data. 

In Fig. |], we plot the relative mass shift be- 
tween the two larger lattices Am = (m(N s = 
16) - m(N s = 20)) /m(N s = 20) for pseudo-scalar 
mesons and the octet baryons as a function of va- 
lence 1 /K va \ for each sea quark mass. We see that 
Am is consistent with zero for the whole range of 
valence quark mass down to the second lightest 
sea quark at K sca < 0.1350. The size effect be- 
comes visible and increases with reduced 1/K va \ 
for the lightest sea quark at K sea — 0.1355. 

The wrapping of valence quarks in the spa- 
tial directions leading to FSE due to the size of 
hadrons is suppressed by the center Z(3) symme- 
try in quenched QCD g. In full QCD, Z(3) sym- 
metry is broken by the wrapping of sea quarks in 
the spatial directions, whose magnitude increases 
toward lighter sea quark. A qualitative under- 
standing of the behavior in Fig. [I] would be that 
the Z(3) breaking turns on rather quickly around 
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Figure 1. Relative mass difference between V = 
16 3 and 20 3 . 

our lightest sea quark. 

The volume dependence of our data is well de- 
scribed by the power law ansatz m(N s ) — m(N s = 
oo ) oc Q. For pseudo-scalar and vector 

mesons, the magnitude of FSE indicated by this 
fit is about 3 % at K va \ — 0.1350, which roughly 
corresponds to the strange quark mass, and in- 
creases to 5 % toward the chiral limit. How- 
ever, since the volume dependence is expected 
to turn into a milder behavior exp[— mps-ZV,] for 
sufficiently large volumes, the actual size of FSE 
should be smaller. In particular, at strange quark 
mass, FSE should not exceed a few % level. 

Finite size effects are more pronounced for 
baryon masses. For our lightest sea quark, the 
magnitude of FSE is about 5 % at K va \ = 0.1350 
and increases to 10% for K va \ = K sca = 0.1355. 
This is comparable with the typical size of the 
quenching error. Therefore, masses of light 
baryons such as ton and toa may be affected by 
FSE. 

3. Hadron masses 

In Fig. g, we compare the valence quark mass 
dependence of the vector meson mass at each 
sea quark mass and in quenched QCD. The full 
QCD data have a steeper slope than the quenched 
QCD, and show a better agreement with the ex- 
perimental points. 

The increase of the slope can be explicitly seen 
in the J parameter ^ shown in Fig. [| The value 
of J, which is consistent with the quenched value 
for heavier sea quark masses, increases as the sea 
quark mass decreases. The chiral limit is still sig- 
nificantly lower than the phenomenological value, 
which may be attributed to quenching of strange 
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Figure 2. Vector meson mass as a function of 
pseudo-scalar meson mass squared. 
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Figure 3. J parameter in full (left panel) and 
quenched QCD (right panel). 

quarks. 

Finite size effects, possibly present at a few per- 
cent level, do not alter these results since masses 
at larger volumes would show a steeper slope in 
Fig.g 

The hadron spectrum at the physical quark 
mass is calculated using data on the largest lat- 
tice. We use and m p to fix the scale and the 
light quark mass m u d- The strange quark mass 
m s is determined from either K or (f> meson mass. 

Our results for the spectrum are plotted in 
Fig. [I| The deviation of the meson spectrum 
from experiment is reduced significantly by the 
inclusion of dynamical up and down quarks. The 
strange baryon masses, such as vn-~ and mo, also 
show a closer agreement with experiment in full 
QCD. However, sea quark effects are less clear in 
lighter baryons like and toa- This is proba- 
bly due to FSE, which is enhanced for the lighter 
baryons. 

We make a comparison of to^-* with if-input 
with the CP-PACS results § in Fig. §. A good 
agreement with their result in the continuum 
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Figure 4. Hadron spectrum with if-input. 

limit suggests that the scaling violation in our 
data is small. 

4. Quark masses 

Calculation of m u d and m s is performed with 
the bare quark mass defined through either the 
vector (VWI) or axial vector Ward identity 
(AWI). For the latter, we use the improved ax- 
ial current with ca at one-loop @. The con- 
tinuum MS quark mass is obtained by the one- 
loop matching R] at scale /i — a^ 1 and evolved to 
/i = 2 GeV using four- loop beta function. 

It is expected that various systematic uncer- 
tainties, such as scaling violation, partially can- 
cel in the ratio defining the AWI mass m q = 
(A 4 P t )/(2(PPt)). This is supported by the good 
agreement with the CP-PACS results in the con- 
tinuum limit in Fig. ||. Such a cancellation 
is not expected in the VWI mass. Indeed there 
is a sizable difference between our AWI and VWI 
results, and scaling violation is large in the CP- 
PACS results. 

We therefore quote the AWI masses in Ta- 
ble | as our central values. The quark masses 
decrease by 20-30% in two-flavor QCD com- 
pared to quenched QCD. Sea quark effects also 
reduce the difference between m s with K- and cj>- 
inputs. These are consistent with the observation 
inRef. || . 

5. Conclusions 

We have demonstrated sea quark effects at 
simulation points in the meson sector using the 
0(a)-improved Wilson quark action and plaque- 
tte gauge action at a -1 ~2 GeV. A closer agree- 
ment of the strange meson and baryon masses 
with experiment is a natural consequence of the 
finding. We also observed that sea quark effects 
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Figure 5. Comparison of wir* (top panel) and 
m s (bottom panel) with CP-PACS results |,|. 
All results are obtained with K-iwput. 



Table 1 

Results of AWI quark masses in MeV units, 
input to fix m s is written in brackets. 
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lead to a significant reduction of quark masses. 
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